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The Ginzburg-Landau model below its critical temperature in a temporally oscillating external 
field is studied both theoretically and numerically. As the frequency or the amplitude of the external 
force is changed, a nonequilibrium phase transition is observed. This transition separates spatially 
1-^ ' uniform, symmetry-restoring oscillations from symmetry-breaking oscillations. Near the transition 

a perturbation theory is developed, and a switching phenomenon is found in the symmetry-broken 
phase. Our results confirm the equivalence of the present transition to that found in Monte Carlo 
simulations of kinetic Ising systems in oscillating fields, demonstrating that the nonequilibrium 
phase transition in both cases belongs to the universality class of the equilibrium Ising model in 
zero field. This conclusion is in agreement with symmetry arguments [G. Grinstein, C. Jayaprakash, 
C/3 ' and Y. He, Phys. Rev. Lett. 55, 2527 (1985)] and recent numerical results [G. Korniss, C. J. White, 

P. A. Rikvold, and M. A. Novotny, Phys. Rev. E (submitted)]. Furthermore, a theoretical result for 
, the structure function of the local magnetization with thermal noise, based on the Ornstein-Zernike 

G I approximation, agrees well with numerical results in one dimension. 
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Bistable systems that are driven between their two states by a periodically oscillating external force are common, 
both in nature and technology. A few examples are hysteresis in ferromagnetic |lj-|^ and ferroelectric ||5|-|7|] materials 
^^ , driven by oscillating applied fields, electrochemical adsorbate systems driven across a phase transition by an oscillating 
^\ . electrode potential [p|-[10[, and liquid crystals driven through a phase transition by pressure oscillations |Q. In this 
^^ ' paper we use magnetic language, henceforth referring to the order parameter as the magnetization and the oscillating 
^^ . force as the magnetic field. 

When the field oscillates at a sufficiently low frequency, the driven system essentially follows the field, switching 

between its two states in a symmetry-restoring oscillation (SRO) with the same period, provided that the amplitude of 

the external force is larger than a critical value which depends on the temperature and the system's spatial dimension. 

S [ At high driving frequencies, on the other hand, the system is unable to follow the field and instead settles down to a 

i-rt ■ symmetry-breaking oscillation (SBO) around one or the other of its zero-field stable states. 

^ ' Over the last decade it has become evident that the boundary between the SRO and SBO regimes corresponds 

O ■ to a singularity that appears to have all the hallmarks of a genuine second-order phase transition. It is therefore 

^ ' appropriate to consider the SRO and SBO as dynamic phases of this far-from-equilibrium system. Characteristic 

J> , features of this nonequilibrium phase transition include a power-law dependence of the amplitude of the SBO on the 

amount by which the frequency fl exceeds its field and temperature dependent critical value flc, as well as critical 

. , slowing-down |12| . The probability density of the period-averaged magnetization exhibits a one-peak structure for 

j^ ' fl < flc and a two-peak structure for Q > Qc p^-p^]. In spatially extended bistable systems, such as the two- 

■ ■ ' dimensional kinetic Ising model below its critical temperature, the transition also displays a divergent correlation 

length and finite-size scaling properties analogous to those familiar from equilibrium phase transitions |l3-n^ . It has 

become common to refer to this symmetry-breaking transition as the "dynamic phase transition" (DPT). 

The DPT was first observed in numerical solutions of a deterministic mean-field equation of motion for a ferromagnet 
in an oscillating field |lq,M, and it has subsequently been seen and studied in numerous Monte Carlo (MC) simulations 
of kinetic Ising systems |[l^- [l7| , p0| -p5| , as well as in further mean-field studies ||l^ , pi| , p3| , P6H27[ . It may also have been 
observed experimentally in ultrathin Co films on Cu(lOO) [0,0. Reviews of earlier research on the DPT and related 
phenomena are found in Refs. |2q,E9|. 



Finite-size scaling analysis of MC data for the DPT in the two-dimensional kinetic Ising model at sub-critical 
temperatures provides strong numerical evidence that this nonequilibrium critical phenomenon belongs to the same 
universality class as the equilibrium phase transition in the two-dimensional Ising model in zero field |l^-[l6[| . While 
this result may seem surprising at first, it is consistent with a symmetry argument by Grinstein, Jayaprakash, and He 
pO| . This argument states that continuous ordering transitions of fully probabilistic cellular automata with Ising-like 
"up-down" symmetry (of which the kinetic Ising model in an oscillating field is an example) should fall in the same 
universality class as the corresponding Ising model in equilibrium. This implies that such a cellular automaton should 
possess an underlying coarse-grained effective Hamiltonian at sufficiently large length scales, which determines its 
universality class. 

The purpose of the present paper is to elucidate the origin of the DPT and to clarify the statistical characteristics of 
the dynamics near the DPT, subject to thermal noise. To this effect we consider a time-dependent Ginzburg-Landau 
model with thermal noise. The equation of motion for the noise-free version of this model is 

^{r,t)^'iJj~i;^ + V^ilj + hcos{nt) , (1) 

where '0(r, t) is the continuous scalar magnetization field, and h and Q are the amplitude and frequency of the spatially 
uniform external magnetic field, respectively. In zero applied field, Eq. ffl) is identical to the conventional Ginzburg- 
Landau equation for the Ising model below its critical temperature. The effects of thermal noise on the system are 
expressed by the stochastic differential equation 

V)(r, i) = V - ^^ + V^V^ + h cos{nt) + R{r, t) , (2) 

where i?(r, t) is a Gaussian white thermal noise. In this paper we show that Eqs. (|l|) and (0) give rise to a bifurcation 
line in the (fJ, h) plane. Near this bifurcation line the equations yield an effective Hamiltonian for a dynamic order 
parameter. This effective Hamiltonian is in the same universality class as the equilibrium Ising model in zero field, 
and its existence provides explicit confirmation of the symmetry argument of Ref. ||30| for this far-from-equilibrium 
system. 

Equations (|l|) and (||) with h — Q give rise to two degenerate ordered solutions only for systems of spatial dimension 
d > 2 at temperatures below criticality. These conditions will be assumed hereafter, unless otherwise is explicitly 
stated. 

The present paper is organized as follows. In Sec. || we show that the spatially uniform oscillation of Eq. (|l|) 
undergoes a bifurcation as O is increased, which separates the symmetry-restoring and symmetry-breaking dynamic 



phases. In Sec. [II we develop a Landau expansion near the bifurcation, which is used to explain the switching 



phenomenon observed in a system subject to thermal noise. In Sec. IV we show that theoretical results for spatial 
power spectra of spin fluctuations (structure functions) obtained by the Landau expansion are in good agreement 
with numerical experiments for a one-dimensional system. A summary and conclusions are given in Sec. M. 

II. BIFURCATION OF THE SYMMETRY-RESTORING OSCILLATION 

In this Section we concentrate on the uniform solutions of the the noise- free system described by Eq. (|y). The 



effects of spatial fluctuations and thermal noise will be discussed in Sec. Ill 
It is easy to see that Eq. (qI) has a spatially uniform oscillation, 

^(t)=^-^^ j^h cos{nt) , (3) 

provided that spatially periodic boundary conditions are used. Without loss of generality, h and 51 are taken as 
positive. Eventually, ip(t) always exhibits a periodic oscillation of frequency fl for any choice of h, f2, and the initial 
value ■0(0)- It exhibits no other periodic or chaotic oscillations. This is so because the dynamical system (^) is 
dissipative and has only two degrees of freedom. 

One should be careful when discussing the dynamics near fl == 0. By shifting time as i ^ t — 7r/(2ri), Eq. (|^) 
reduces to 

TJ;(t) = ^ - tP^ + h sin(m) . (4) 

If one puts f7 = in Eqs. (0) and (0) while keeping h finite, they have different fixed points. The period T (= 2Tr/fl) 
of the applied field tends to infinity as fi — *■ 0. One should therefore discuss the long-time behavior oi ip {t ^ T) at 
finite fl, and then take the limit fi — > 0. The above discrepancy originates from the interchange of the limits t -^ oo 
and ri ^ 0. If one takes the limits correctly, the long-time behaviors of Eqs. (|^) and (Q) give the same results. 



For h = 0, tfj{t) eventually approaches one of the stable fixed points, ^o = il- Then, under an applied field 
hcos{p,t) with a small amplitude, it is easy to see that ip{t) exhibits a periodic oscillation. In fact, to first order in h, 
Eq. @ is solved by 

V'(i) = ±n -2 [2 cos{nt) + n sm{nt)] (5) 

for t -^ oo. This is a SBO since TJj{t) oscillates near either tp ~ +1 or —1, depending on the initial condition. We thus 
expect that Eq. (0) exhibits a symmetry-breaking periodic oscillation in the regime of relatively weak h. 
Let V'CO be a solution of Eq. (^). It is then easy to show that ■0(t) given by 

m = -v- (t + 1) (6) 

is also always a solution of Eq. (g), including even the transient process. As a special case (see Eq. ( p7| ) below), Eq. (^) 
has a particular solution with the symmetry 

V'(i)=-^(i+|). (7) 

If the dynamical behavior satisfies the symmetry (m) , one obtains 

i-T 

i}{t)e^^^^dt = , (^ = 0, ±2, ±4, • • •) . (8) 



However, the fact that the system (0) has the symmetry (M) does not necessarily mean that the dynamical behavior 
always exhibits this symmetry. In fact, as discussed above for small h and as shown below, Eq. (j^) has a stable 
symmetry-breaking solution for a certain range of h and f2. 

As discussed above, the dynamics in a weak external field shows a SBO. This implies that the SRO, if it exists, 
should do so for a relatively large amplitude of the external field. This also suggests that there should exist a 
transition between the SBO and the SRO, provided that the SRO stably exists. This is an immediate consequence of 
the symmetry consideration. 

For the moment, let us consider the parameter values h = 1.0, and fi = 1.08 and 1.1. For these parameter values, 
the system has attractors as shown in Fig. nl Throughout this paper, the numerical integration of Eq. (0) is carried 
out by using the fourth-order Runge-Kutta algorithm with the time increment At = T/1024 for all frequencies. The 
attractors are limit cycles of period T. They have the symmetry (0) for fi = 1.08, but are asymmetric for Q = 1.1. 
The above considerations suggest the existence of a phase transition between these different characteristic oscillations. 
Figure g shows the hysteresis loops, i.e., the dependence of ^{t) on h{t) = hcos{flt) for Q below and above flc, the 
critical frequency separating the interchange of the symmetric and non-symmetric oscillations. Numerically, we find 
ric w 1-095 for h = 1.0. 

Next we consider the stability of the attractor with the symmetry (R), shown in Fig. 1^, as ft is increased at fixed 
h. The stability of a periodic oscillation is discussed as follows. Let ^lj{t) be a particular solution of Eq. (||) on an 
attractor, which may be either stable or unstable. In order to examine its linear stability, we seek the temporal 
evolution of the deviation S^{t) from this solution. Then, 6^p{t) obeys the equation of motion, 

Sm - [1 - 3iij{t))^]S^{t) . (9) 

Since 'ip{t) is periodic with period T, Sijj{t) is solved as 

Sipit) = B{t)e^'Si;{0) . (10) 

If we define 

A = l-3V^=l-3^y msjfds, (11) 



then 

B(t) = exp 



-3 A(^(t'))'-^'K 
Jo 



(12) 



is a periodic function of period T, i.e., B{t + T) = B{t). Here we have defined the period-average of f{t) as 
7(i) = T-i J^ f(t + s)ds. The results (|lO| - ^ follow from the Floquet theorem |3^]. The quantity A is called the 



Floquet exponent and indicates the stability of the periodic oscillation under consideration, i.e., ip{t) is linearly stable 
(unstable) if A < (> 0). Numerical results for A calculated by Eq. (0) are shown in Fig. 0. For fi below Oc, the 
critical value for a given h, A takes a negative value, which is denoted by A. The limit cycle for 17 < fie is symmetric 
as shown in Figs. Ma) and 0(a). As fi is gradually increased, the Floquet exponent approaches zero and again takes 
a negative value for ^ > flc. Figures |l|(b) and ||(b) show the stable attractors (solid curves) corresponding to the 
limit cycles for 51 > flc- Figure H shows the stability regions of the SRO and the SBO. In the SBO region, one finds 
that there exist two attractors, C- and C+, one of which is chosen depending on the initial condition. For fl > flc, 
there is also a symmetry-restoring unstable limit cycle, whose Floquet exponent is denoted by A„ in Fig. 0, and whose 
trajectory is depicted by the dashed curves in Figs. |](b) and||(b). The transition at Qc is continuous, as is expected 
from the frequency dependence of the Floquet exponent shown in Fig. y. 

The unstable limit cycle, i.e., the SRO for O > flc, is numerically obtained as follows. Taking an initial value ip„ 
at time i„ = nT, then integrating Eq. (|3|) until tn+i = tn + T, we obtain tpn+i- In this way we get the ipn+i vs ipn 



curve, 



■)pn+l = gi'ipn) ■ (13) 

Examples of numerically obtained g{ij) are shown in Fig. H. Figure 0(a) is for Q < flc, and Fig. ra(b) is for Q > flc- 
Depending on fi, there are one or three fixed points tpf satisfying ■0/ = dii^f), which correspond to cross-sections of 
limit-cycle attractors. The stability of a limit cycle is determined by the slope of g{ip) at ip = ipf, i.e., the Floquet 
exponent is given by 

A=lln|5'(^/)|. (14) 

The unstable periodic orbit shown in Figs. |lj(b) and y(b) is the one numerically integrated with the initial value i/jm, 
the unstable fixed point. The temporal evolutions of one unstable and two stable oscillations arc shown in Fig. ra. 
One should note that if the stroboscopic map is constructed for times i„ = t + nT, the form of g{4>) depends on r. 
However, the number of fixed points of i/'n+i = gi^^n) and the corresponding slopes, which yield the Floquet exponents 
for the fixed points, are independent of r. 

The bifurcation point flc depends on h. The theoretical bifurcation curve given by the solid curve in Fig. is 
determined as follows. We first expand ■ip{t) as the Fourier series 

oo 

V^(t) = Y. Mt)e''''' , (15) 

£=-00 

where ip-i — ip^. The temporal evolution of the coefficients {ipi{t)} is assumed to be much slower than the time scale 
T. Inserting Eq. dl3) into Eq. {m and comparing the coefficients on both sides of the equation, we obtain 

ipe + imiJi = V(? - X! X! i'm^n^e-TJi-n + 2 (^^,1 + ^e-i)h ■ (16) 

m n 

From the symmetry argument, Eq. (m may have a solution with the symmetry (M). If the limit cycle under consider- 
ation is symmetric, we find from Eq. (j^ that 

i^i^O for£ = 0,±2,±4,--- . (17) 

We now consider the stability of this symmetric oscillation. As the simplest non-trivial approximation, we use the 
truncation i = and ±1, which yields 

i^o=[l-6|7/^ip-Vo']^o, (18) 

V^i + zflVi = [1 - 3|7/^iP - 3^g] ^i + ^h. (19) 

The above equations have a SRO (■00 — 0) provided that 1 — 6|V'iP < 0. On the other hand, for 1 — 61-01 P > 
the steady-state value of -00 does not vanish, which implies the emergence of a SBO. Therefore, we find that the 



boundary between the regions of stability of the SRO and SBO is determined by IV'^P — ^- Combining this with the 
steady-state value of ipT obtained from Eq. (jlj) , the bifurcation point fl^ for fixed h is determined by 



Qr 




or, equivalently, 




(20) 



(21) 



One finds that this kind of bifurcation is observed for h larger than a critical value, l/v6 in the above approximation. 
The curve given by Eq. (21) corresponds to the transition line, which in Fig. His compared with results from numerical 
integration of Eq. (|^). For il < fie, there is only one type of periodic motion, namely the symmetry-restoring one. For 
il > ric, on the other hand, there are two types of oscillations: one is the unstable SRO, and the other is represented 
by the two stable SBOs, one of which is observed for a given initial condition. 

Here a comment on the critical value of h under a static field, f2 = 0, should be added. The above approximation 
yields the critical value l/-\/6 « 0.408. On the other hand, the standard calculation in mean-field theory yields the 
spinodal field as the field where the metastable minimum in the ip^ potential disappears. This condition requires that 
the equations. 






(22) 



are simultaneously satisfied. The second equation gives V'spinodai — ±l/\/3 which, when inserted into the first equation, 
yields /igpinodai = 2\/3/9 ~ 0.385. This value is about 6% below that obtained in the above discussion. We carried out 
numerical calculations for values of 51 as small as 0.05. The numerical results seem to be closer to /ispinodai than to 
the present approximate value, l/\/6. However, calculations at even smaller fl, which were not feasible in the present 
study, would be needed to reach a firm conclusion. A sharp decrease of the critical value of /i as $7 is decreased may 
suggest the possibility that the transition curve may have a kind of singularity, i.e., that dh/dfl might diverge as f2 
approaches zero p^ . 

Next we evaluate how the amplitude of the SBO develops for O above ilc- The steady-state values V'o'^ and ipf^ are 
obtained by setting ipQ^ — sin 9 and V'r = 6~2 cos 9 ■ e*". For fl < il,c, 9 — 0, while 9 is small for $7 > fie- A short 
calculation shows that the order parameter ipQ^ is asymptotically given by 



V-o'^iciV^-rfc 



with ci = ^/MJiAn 



2 -L 11). The amplitude and phase of ipf are given as 

1 



iV'r 



1 



i~-ci{n~nc 



and 



respectively, where 



a = ac + C2{i^ ~ ^c) 



(23) 



(24) 



(25) 



smttc 



C2 = 



-1 


V4^? + 1 ' 


22 



4172 + 11 



(26) 



The 51 dependences of the amplitudes of the SRO and SBO are shown in Fig. M, which was obtained from the stable 
and unstable fixed points of Eq. (nSl) . 



Numerical integration shows that the time evolutions of the symmetry-breaking orbits 0+ and C_, which are ip+{t) 
and ip- (t) , respectively, are related as 

^p.(t)^-xb+{t+^) (27) 

(see Fig. ^), where 

V-i (t) = V'i (t) - (V'i (t))^ + h cos{nt) . (28) 

The symmetry ( p7| ) is just a particular case of the general symmetry relation (o). 

We find that the stable limit cycle of Eq. (||) for fi < Oc (SI > Sic) is a symmetry-restoring (symmetry-breaking) 
oscillation, and that one of the two stable SBOs is chosen, depending on the initial condition. Next, we examine 
the stability of the uniform (stable) oscillation with negative Floquet exponent against inhomogeneous fluctuations. 
(Remember that the solution of Eq. (pf) is the uniform solution of Eq. (n|).) Let tp{t) be a stable solution of Eq. (H), 
which implies that its Floquet exponent is negative, i.e., SRO for SI < Sic or SBO for SI > Sic. The negative Floquet 
exponent is denoted by A instead of A. Next, let ip{r,t) be the deviation from ijj{t), i.e., 

V'(r,i)=^(i)+V^(r,i). (29) 

The Fourier transform of 7/'(r, t) obeys 

Mt) = [1 - 3(VXt))'] M-t) - k'Mi) , (30) 

provided that the deviation is sufficiently small. This equation is solved as 

^k(i) = Bit)e^^*MO) , (31) 

where B(t) is again a periodic function and 

Ak = A - fc^ (32) 

is the linear growth rate of the Fourier mode at wave vector k. Since A < 0, Ak is always negative, which implies that 
the uniform oscillation with negative Floquet exponent is linearly stable against inhomogeneous fluctuations with any 
wave vector. This implies that the system (^J) eventually approaches a spatially uniform oscillatory motion, provided 
that there exists no other stable dynamical behavior. 

For simplicity, the values of Qc{h), obtained above from the spatially uniform solution, will be referred to as the 
mean-field values of Sic. In spatially extended systems with thermal noise, the actual values of Sic are renormalized 
by fluctuations. 

III. LANDAU EXPANSION AND THERMAL NOISE EFFECTS 

We now move on to the discussion of the spatially extended system with local interactions, which is described by 



Eqs. M) and (Eh. The noise- fr ee cas e, Eq. (ry), is discussed in Sec. Ill A. The effects of thermal noise, described by 



Eq. (g), are considered in Sec. |[IIB 



A. Landau expansion near the bifurcation point 

Let ilj*(t) be the spatially uniform SRO which obeys Eq. (0) and satisfies the symmetry (1^), and let A be its Floquet 
exponent, which is calculated by Eq. (O) with ^{t) — 4'*it)- It is given by A = A(< 0) for SI < Sic, and K = \u{> 0) 
for SI > Sic, in the notation used in FigTH. Expanding ■(/'(r, t) around this SRO as 

i!{Y,t)=i,,{t)+BS)(p{Y,t) , (33) 

where i3*(i) is defined by Eq. (|2|) with ijj{t) = V'*(^), and inserting this into Eq. (|]), we immediately find 



(r,<) = (A + V^ 



2j,3 



37/.*(i)S,(i)0^-(S,(t))^0' 



(34) 



Note that because of the particular symmetry (0) in the SRO phase, _B*(t + T/2) = B^,{t). Since the coefficients of 
the above equation are periodic in time, we may use their time-averaged values, noting that the characteristic time 
|A|~^ of near the transition is much longer than T. Making use of the symmetry relation ([T]), one can prove that 



- / Mt)B*{t)dt=- / V*Wexp 



Thus, Eq. 



T 



reduces to 



T 



(r,t) = (A + V2 



-3 / {{M'S)f - ^l}ds 





^3 ^ ^ SH{c^} 



dt = 



(35) 



(36) 



with b = Bl and 

Here, F is a positive constant, which will be determined in an appropriate way below. 



±^^ + \iy^f + \^^ 



dr 



(37) 



B. Switching phenomenon 



The thermal noise effects near the DPT were studied in Refs. O-IT^. In the present continuum model, the thermal 
noise i?(r, t) is included in Eq. (0) as a Gaussian white noise satisfying 



(i?(r,f))=0 
{R{r, t)R{v', t')) ^ 2V5{r - v')5{t - t') 



(38) 



where (• ■ ■) denotes the ensemble average. The noise strength T is proportional to the temperature of the system. 
In Refs. |13[-p6[ the thermal noise effects were studied by observing the time evolution of the total magnetization 
in two-dimensional kinetic Ising systems. A switching phenomenon between two asymmetric oscillatory states was 
observed for values of f2 slightly above fJc- The origin of this phenomenon in the present continuous-spin model (0) 
can be understood as follows. Inserting the expansion (p3) into Eq. (0), and approximating the coefficients by their 
time averages, we obtain 



(r,t) = (A + v2 



+ /(r,t) 






fir,t), 



(39) 



where /(r,t) = (-B*(i)) ^R{r,t) is a Gaussian white noise with the strength T = B^ ^T. This is also chosen as the 
value of f in Eqs. (^) and (|3). 

Equation (p3) is identical to the conventional c/)^ Ginzburg-Landau equation in zero external field with a thermal 
noise term. This equation belongs to the same universality class as the Ising model [Q. This is the mechanism of 
the DPT and the switching phenomenon observed in Refs. ]13-|lffl. 

Equation ( |39| ) is the central result of this paper, which makes the connection to previous work on the DPT in 
kinetic Ising models. In those studies, the local dynamic order parameter has been taken as the period-averaged 
magnetization. 



Q„(r) ^■iP{r,t„) = - 



tn+T 



ip{r,t)dt 



(40) 



where n — mt(t/T). Since (/)(r, t) depends on t only on time scales much longer than T, it can be replaced by a variable 
0„(r). Thus it is easy to show that the traditional form of the local dynamic order parameter is simply proportional 
to (j)n{r): 



Qn{r) = ^ 



tn+T 



t-n 



B{s)(t){Y,s)dsK. — 



B rtn+T 



{r,s)ds = B(t>n{r) . 



(41) 



w 



The global dynamic order parameter is simply the spatial average of Qr^ 
(/)(r, t) are also proven for the traditional dynamic order parameter, (5„. 
The Fokker-Planck equation corresponding to Eq. (^9|) takes the form 



(r). Thus, any results that are proven for 



[ s r s f p{cj^,t} \ 

J <50(r) [ ^^Ucj^ir) \P*mJ 



dv 



where 



P*{(t)} oc e-^'-f-^J 



(42) 



(43) 



is the steady-state probability density, which has a single (double) peak structure for 17 < fie (^ > ^c)- Here 7i*{0} is 
the single- or double-peaked renormalized potential function. For fi > ri^i the well separation in 7i*{0} is proportional 
to (yi — ^c)^ , where j3 is the magnetization exponent for the Ising model in the appropriate spatial dimension [0-n6|. 
We postulate that the dynamics of the total magnetization (per unit volume, V)^ 



Mt) 



1 

V 



{r,t)dr 



(44) 



for Eq. (p9p, takes approximately two values. This implies that the dynamics can be modeled by the Langevin equation 

Mt) = A0o(t) - b{Mt)f + fo{t) , (45) 

where fo{t) is a random force with 

ifoit)) = 
{foit)foit'))^2ToSit~t'). (46) 

The Fokker-Planck equation is thus approximated by 

. d 



l^(^-^)-^«4 



p,. , . ^ f P{cbo,t) 



where 



io) oc 



-■HUM 



where 7io(0o) corresponds to the critical order-parameter distribution for the Ising model 
absence in the volume-averaged 0o of spatial variations, this result is analogous to Eq 

As an illustration of the switching behavior, Fig. g shows the evolution of the tota 
numerically solving Eq. (ph for d = 1 with fi slightly larger than the mean-field value of fl 
of Eq. (||) was carried out by the second-order stochastic Runge-Kutta (SRK) algorithm 
lattice points with lattice spacing Ax (set to 0.5 throughout this paper) and using 



(47) 



(48) 

5|,||,||. Except for the 
_ for H* ((/)). 

magnetization obtained by 
. The numerical integration 
B3] , dividing the space into 



(V^^), = 



(Ax)2 



(49) 



at the lattice site j with periodic boundary conditions. A switching phenomenon is clearly observed. 

The switching phenomenon can be formulated in a different way as follows. The temporal evolutions '0+(<) and 
V'-(i) of the symmetry-breaking orbits, C-|_ and C_, respectively, obey Eq. (|2^). We study the additive noise effect on 
the dynamics, adding a weak noise R, as in Eq. (0). We also include the spatial variation of the dynamical variable, 
adding V'^ip. Consider a local bistable variable. For fi > ilc, depending on the initial condition, either C+ or C_ 
is selected, provided the noise is absent. If the noise is sufficiently weak, the phase point is almost always on either 
C+ and C-. When V'+ and ip- are close, the phase point can switch to the other orbit through the noise effect. The 
above picture can be mathematically formulated as follows. Let a{t) be a variable which takes the values -1-1 (—1), 
provided the phase point is on C+ (C_) at time t. Then the temporal evolution of a{t) is approximately described 
by the switching dynamics. The above picture is generalized by introducing the position dependent variable a(r,t), 
defined by 



V'(r,i) = 



l-a(r,i) 



^+it + 4) 



l + a{r,t) 



i^-it + to) 



(50) 



where t^ are certain initial times. The variable a(r, t) indicates whether the local magnetization ■!/)(r, t) is close to 
V'+(i + iJ) or V'-(i + to')- Namely, if a(r, i) is near +1 or —1, ip{r,t) is close to V'+(i + io') or ?/'-(i + io'), respectively. 
Without loss of generality, ig are chosen such that cos{ntQ ) = 1, and therefore we put t^ = 0. Inserting Eq. ( [50| ) 
into Eq. (0), after some algeSra we rigorously get 



1 



\2 „ off. I. /i\\2 / /. f,\\2 



d(r,t)--(l-a^)^(^+(t)-V-(t)) a-3{(V'+(t))'-(V-(t))nj+V^a + g(r,i) (51) 

with 

, , -2R(r,t) 

Furthermore, as long as the thermal noise is weak, the average switching time between a — +1 and —1 is long. The 
temporally periodic coefficients can therefore be replaced by their average values, which reduces Eq. ( |5l|) to 

d = ^(l-a^)a + V^a + g(r,i) (53) 

with 

M = J(V'+ W - ^-(t)? = l{B{t)y (0+(t) - 0_(i))' (> 0) , (54) 

where -0^ = ^^ by symmetry. The quantities B{t) and (t)±{t) are the same as in Sec. |l[ Equation ( p3| ) again takes 
the form of the Ginzburg-Landau equation with a double- well potential, with stable fixed points a = ±1, provided 
that the spatial variation of a and the noise are neglected. Equation (p3) shows the switching phenomenon. 

IV. STRUCTURE FUNCTION IN THE SRO PHASE 

In this Section, we study the structure function in the SRO phase (A < 0) in a one-dimensional system. We define 
the structure function, 



Ik = {\Mt)n , (55) 



for the Fourier transform ^Jk{t) of the period-averaged order parameter. 



I rt+T 



^lT,t) = -j ,p{r,s)ds. (56) 

We obtain /k by numerically solving Eq. (||) in a one-dimensional system for several values of Q in the SRO phase 
as shown in Fig. H. One characteristic feature of the numerically obtained structure function is the k~ behavior 
observed for relatively large k. Here, Ik can be evaluated using the Landau expansion with thermal noise as follows. 
By inserting Eq. (BSh into Eq. (p^) and using the approximation [see also Eqs. ( pO| ) and (|4T 

I i-t+T -g rt+T 



i^ir, t)^-j Bis)<Pir,s)ds^-j q^{r,s)ds, (57) 



the structure function is evaluated as 



I^. = ^j dsj ds'{Ms)K{s')). (58) 



Using the linearized form of the equation of motion for (/)(r, t), Eq. (B9f), which yields 

</)(r,i) = (A + v2)</,(r,i)+/(r,t), (59) 

we immediately get the correlation function 

{Ms)Kis')) = -e-^-^l^-^'l , (60) 

7k 



where 

u2 



7k = |A|+A;^ (61) 

is the dampmg rate for the fluctuation 0k (t) ■ Substitution of Eq. ( |60| ) into Eq. ( p8| ) yields 

/k = prrG(7kr) , (62) 



where p = B B '^ and 



2 / 1 — p~^' 

G(^) = z? 1-^:— ) (63) 



is a scaling function. We thus find that the structure function /k behaves asymptotically as Jk ^ 7k for 7kT <^ 1 
and as ^ 7k ^-^""'^ ^^^ 7kr ^ 1- In the former case, the relation Ik (x 7^^ is identical to the Ornstein-Uhlenbeck form 
which holds for T ^ 0. The latter characteristic explains the k~^ behavior in the large-fc regime. This unexpected 
behavior originates from the temporal averaging procedure (pq), which ensures that the interfaces between regions of 
positive and negative values of V'k(i) are not thin. As a result, /k does not obey "Porod's law," /k ~ k~'^'^^^^ [ pq , 
which would yield /k ~ fc~^ for the present case of d = 1. 

The solid curves in Fig. O correspond to the theoretical result (p2-63). Although this linear model for the fluctuations 



agrees quite well with the numerical simulations, particularly in the large-wavenumber regime and for strong thermal 
noise, nonlinear fluctuations play a significant role and the linear model eventually breaks down. 

V. SUMMARY AND CONCLUSION 

In this paper we used a time-dependent Ginzburg Landau model in a temporally oscillating external field to 
understand the dynamic phase transition (DPT) observed in Monte Carlo simulations of the corresponding kinetic 
Ising model below its critical temperature [l^- p7| , |2l| - p^ , p8| , p9[ . Analyzing the stability of spatially uniform oscillations, 
we found a bifurcation of the symmetry- restoring oscillation (SRO), which leads to the onset of a symmetry-breaking 
oscillation (SBO). Developing a Landau expansion near the bifurcation point, including additive thermal noise, we 
obtained an effective Ginzburg Landau Hamiltonian for the amplitude of the SBO, which is proportional to the 
dynamic order parameter which characterizes the DPT. This effective Hamiltonian has the same form as the standard 
(/)■* Ginzburg Landau Hamiltonian in zero external field, which describes the long-range properties of the Ising model in 
the critical region p^. This result implies that the DPT belongs to the same universality class as the equilibrium Ising 
model in zero external field, in agreement with recent high-precision numerical results from Monte Carlo simulations 
p3|-p6[, as well as with a symmetry argument which states that the equilibrium Ising universality class should 
encompass all stochastic cellular automata with Ising "up-down" symmetry [ pO[ . 

To the best of our knowledge, the work presented here is the first in which an effective Hamiltonian has been 
explicitly derived for a far-from equilibrium phase transition, confirming that the transition belongs to the same 
universality class as the equilibrium Ising model in zero field. The result represents a significant expansion of the 
realm of validity of symmetry arguments from equilibrium to nonequilibrium phase transitions. 
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Fig. 1 



FIG. 1. Limit-cycle attractors for parameter values (a) h — 1.0, Q — 1.08 and (b) h — 1.0, Q — 1.1. The phase points move 
clockwise. In (a) {Q < Qc ~ 1.095), there stably exists only one limit cycle, which is symmetric in the sense that Eq. (pi) is 
satisfied. In (b) (SI > Qc), the symmetric limit cycle denoted by the dashed curve is unstable, and there appear two stable 
non-symmetric limit cycles, C+ and C_. 
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Fig. 2 



FIG. 2. Hysteresis loops of limit-cycle attractors, i.e., tp{t) vs h{t) = hcos{Q,t) for (a) h = 1.0, Q = 1.08 and (b) h = 1.0, 
SI = 1.2. The phase points move counterclockwise. In (a) (f2 < fie « 1.095), there is only one stable symmetric limit cycle. In 
(b) {Q > flc), the symmetric limit cycle (dashed curve) is unstable, and there exist two stable non-symmetric limit cycles, C+ 
and C_. 
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Fig.3 

FIG. 3. The Floquet exponent A, shown vs 57 for h = 1.0. For Q. < Qc, the Umit cycle is symmetric and stable. For Q > Jlc, 
the dashed line is the Floquet exponent for the unstable symmetric limit cycle, and the sohd line is the exponent for the stable 
symmetry-breaking limit cycle. For details on the calculation of Floquet exponents, see the text. 
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FIG. 4. The bifurcation curve separating the symmetry-restoring oscillation (SRO) and the symmetry-breaking oscillation 
(SBO). The numerically obtained points are represented as data points, and the approximate theoretical result (gl|) as a solid 
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Fig. 5 



FIG. 5. Stroboscopic maps tp„ = ip{n ■ 2-k/Q.) for (a) Q. = 1.08 (< Q.c) and (b) Vi = 1.2 (> Q.c) with h = 1.0. For VI > 0.^ 
there appear two stable, non-syminetric limit cycles. 
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FIG. 6. Temporal evolutions of the stable orbits, C+ and C- (solid curves), and the unstable orbit (dashed curve) in the 
SBO regime for h = 1.0 (< Q.c) and Q. = 1.2 (> fie). The dot-dashed curve represents h{t). 
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FIG. 7. Bifurcation diagram of ipn{= i^in ■ Stt/H)), i.e., the fixed points of the map ( |13[ ) for h = 1.0. For SI < fie ~ 1.095, 
there is one unique fixed point. For Q > Qc, there exist one unstable fixed point (dashed line) and two stable fixed points (solid 
curves) . 
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FIG. 8. Switching phenomenon generated by Eq. (H) in a one-dimensional system for Q slightly larger than the 



-{16T) 



mean-field value of flc- Here M{t) — L ^ ( ip[x,t) dx is the total magnetization, where L is the system size and 



-(mT) 



ft + niT 



ip{x,ty ' — (mT)"^ J^'""^ ■4'i^jS)ds. Parameters are h = 1.0, fl = 1.2, [Qc ~ 1.095) and F — 0.005. System sizes are 
(a) L = 16Aa;, (b) L — 20Aa;, where Ax = 0.5 is the lattice spacing. The time increment is chosen as At = r/1024 = 0.005113. 
Numerical simulations were carried out for L/Ax — 64, 80, 96, 112, 128, 144, 160, 176, and 192. The average time between 
switching events was observed to increase monotonically with L. 
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FIG. 9. The structure functions for the Fourier components of '(/'(r, t) in a one-dimensional system with thermal noise in the 
SRO regime. The numerically obtained results are shown as o, while the theoretical result of the linearized model (M) is shown 
as a solid curve. The insets are log-log plots. The parameters are (a) system size L — 8192 x Ax with lattice spacing Ax — 0.5, 
h = 1.0, D. = 1.0, (A = -0.220), and T = 0.05, and (b) L, Ax, h, and F the same as in (a), while Q. = 1.08, (A = -0.0307). 
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